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Abstract : We investigate the limit behaviour of the spectral measures of matrices following the Gibbs 
measure for the Ising model on random graphs, Potts model on random graphs, matrices coupled in a chain 
model or induced QCD model. For most of these models, we prove that the spectral measures converge 
almost surely and describe their limit via solutions to an Euler equation for isentropic flow with negative 
pressure p{p) = — 
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1 Introduction 

It appears since the work of’t Hooft that matrix integrals can be seen, via Feynman diagrams expansion, as 
generating functions for enumerating maps (or triangulated surfaces). We refer here to the very nice survey 
of A. Zvonkin’s . One matrix integrals are used to enumerate maps with a given genus and given vertices 
degrees distribution whereas several matrices integrals can be used to consider the case where the vertices 
can additionally be coloured (i.e can take different states). 

Matrix integrals are usually of the following form 

Zj,iP) = I 

with some polynomial function P of d-non-commutative variables and the Lebesgue measure dA on some 
well chosen ensemble oi N x N matrices such as the set Hn (resp. Sn, resp. SympN) oi N x N Hermitian 
(resp. symmetric, resp. symplectic) matrices. One would like to understand the full expansion of Zi\f{P) 
in powers of N. For instance, in the case where the matrices live on TdjV) the formal expansion linked with 
Feynamn diagrams is of the type 


^,o,Z,{P) = J2^CAs) 

_ 9>0 
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where Cp{g) enumerates some maps with genus g. Such an expansion was proved to hold rigorously in the 
one matrix case by K. McLaughlin and N. Ercolani in 2002. 

A related issue is to understand the asymptotic behaviour of the corresponding Gibbs measure 


g^{dA^ ■ ■ ■ dA^) = ,A^a))dA^ ... dA 


Zn{P) 


N 

d * 


More precisely, iffora A^xA^ matrix A, (Ai(A), 


, Xn^A)) denotes its eigenvalues and := A^ ^ Si=i 


its spectral measure, one would like to understand the asymptotic behaviour of , • ■ • , ) under the 

Gibbs measure g.p when N goes to infinity. Of course, this understanding is intimately related with the 
first order asymptotic of the free energy Fn{P) = N~‘^ XogZ^iP)- In fact, the rigorous approach of the full 
expansion of matrix integrals when d = 1 given by K. McLaughlin and N. Ercolani is based on Riemann 
Hilbert problems techniques which themselves require a precise understanding of such asymptotics of the 
spectral measures. 

However, only very few matrix integrals could be evaluated in the physics litterature, even on a non 
rigorous ground. These cases corresponds in general to the case where integration holds over Hermitian 
matrices. Using orthogonal polynomial methods, Mehta obtained the limiting free energy for the Ising 
model on random graphs, corresponding to d = 1 and P{A, B) = P{A) + Q{B) — AB when P{x) = Q{x) = 
gx'^ + . He extended this work with coauthors to matrices coupled in a chain, model corresponding 

to P{Ai, ■ ■ ■ ,Ad) = Pi{^i) — However, he did not discuss in these works the limiting 

spectral distribution of the matrices under the corresponding Gibbs measure. On a less rigorous ground, P. 
Zinn Justin |^, discussed the limiting spectral measures of the matrices following the Gibbs measure 
of the so-called Potts model on random graphs, described by P{Ai, • • • , Ad) = Pii^i) ~ J2i=2 
Very interesting work was also achieved by V. Kazakov (in particular for the so-called ABAB interaction 
case), A. Migdal and B. Eynard for instance. We refer to the review |Q of B. Eynard for a general survey. 
Matytsin obtained the first order asymptotics for spherical integrals, from which he could study the 
phase transition of diverse matrix models (see iQ for instance). O. Zeitouni and myself |]^ gave a complete 
proof of part of his derivation in fl^ and the present paper is actually finishing to put his article |18 on a 
firm ground. 

In this paper, we investigate the problem of the first order asymptotics of matrix integrals with AB 
interaction, including the above Ising model, Potts model, matrix model coupled in a chain and induced 
QGD models. The integration will hold over either Hermitian matrices or symmetric matrices. The case of 
symplectic matrices could be handle similarly. We obtain, as a consequence of [T^ , the convergence of the 
free energy and represent its limit as the solution of a variational problem. We here study this variational 
problem and characterize its critical points. One of the main outcome of this study is to show that under 
the Gibbs measure gp of the Ising model described by 


P{A, B) = P{A) + Q{B) - AB 

with P{x) > ax'^ -I- b and Q{x) > cx^ -I- d with a, 5 > 0, the spectral measures of (A(^, ) converges almost 

surely and to characterize its limit. More precisely, we shall prove that 
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Theorem 1.1 1) converges almost surely towards a unique couple (hajII'b) of probability measures 

on ]R. 

2) {p,A, P-b) are compactly supported with finite non-commutative entropy 



y\dp{x)dp{y). 


3) There exists a couple {p 


A^B A- 


of measurable functions on IR x (0,1) such that {x)dx is 


a probability measure on IR for all t € (0,1) and {pA, Pb, P^^^are characterized uniquely as the 


minimizer of a strictly convex function under a linear constraint (see Theorem 3.S). 

In particular, {p^^^, are solution of the Euler equation for isentropic flow with negative pressure 


p{p) = —^p^ such that, for all {x,f) in the interior of Tl = {{x,f) € IRx [0, l]',p^^{x) ^ 0}, 


\ dt{pt-^ut-^) + - ^{pt-^f ) = 0 

with the probability measure pf^^{x)dx weakly converging towards pA{dx) (resp. psidx)) as t goes to zero 
(resp. one). 

Moreover, we have 

P'{x) - X - ^Uq^^{x) - ^IIpa{x) = 0, and Q'(x) - x- ^Hpb{x) = 0 
in the sense of distributions that for all h G Cl{IR), 

f J dpA(,x)dpA{y) = J {P'{x) -X - ^u^^^{x))h{x)dpA{x), 

7 [ - ^^^dpB{x)dpB{y) = [iQ'ix) - X + ^u^^^{x))h{x)dpB{x). 

4: J x-y J 2 

A more detailed characterization of (pA, Pb, P^^^is given in Theorem \3.t\. 


Here, Hp stands for the Hilbert transform of the probability measure p given by 

(x-y) . . 

jdp{y) 


IIp(x) = PV [ —-— dp{y) = lim [ ■ 
J x-y eio J 


(x-j/)2 , 

To obtain such a result, we shall first study the limit obtained in ||^ for spherical integrals. This limit was 
indeed given by the infimum of a rate function over measure-valued processes with given initial and terminal 
data. We show in section ^ that this infimum is in fact taken at a unique probability measure-valued path, 
solution of the Euler equation for isentropic flow described in (pTI]). Using a saddle point method, we derive 
from in Theorem 3.1 formulae for the limiting free energy of some matrix models with AB interaction. 
In the Ising model case, this free energy is indeed written as the infimum of a strictly convex function, from 
which uniqueness of the minimizers is obtained. As a consequence, we obtain the convergence of the spectral 
measures under the Gibbs measure for Ising model. A variational study then shows that the limiting spectral 
measures satisfies the above set of equations (see Theorem |3.3|) . For the other considered models (q-Potts 
model, matrix coupled in a chain, induced QCD), obvious convexity arguments and therefore uniqueness is 
lost in general, but still holds in certain cases. However, we can still specify some properties of the limit 
points (see Theorem 3.5). 
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In this paper, we shall denote C([0, 1\,V{IR)) the set of continuous processes with values in the set V{Si) 
of probability measures on M, endowed with its usual weak topology. For a measurable set of IR x [0,1], 
denotes the set of real-valued functions on which are p times continuously differentiable with 
respect to the (first) space variable and q times continuously differentiable with respect to the (second) time 
variable with bounded derivatives. C^’'^(n) will denote the functions of with compact support in 

the interior of the measurable set fl. Lp{dii) will denote the space of measurable functions with finite 
moment under a given measure fi. We shall say that an equality holds in the sense of distribution on a 
measurable set fl if it holds, once integrated with respect to any functions. 


2 Study of the rate function governing the asymptotic behaviour 
of spherical integrals 


In 1^, Ofer Zeitouni and I studied the so-called spherical integral 

if {Dm, Em) := ^ exp{Vtr(C/i?wC/*£;^)}dm^(t/), 

where irf denotes the Haar measure on the orthogonal group Om when /3 = 1 and on the unitary group 
Um when /3 = 2, and Dm, Em are diagonal real matrices whose spectral measures converge to We 

proved (see Theorem I.l in [|^) the existence and represent as solution to a variational problem the limit 

lim N-Hoglf{DM,EM). 

N^OO 

This result in fact was obtained under the additionnal technical assumptions that there exists a compact 
subset K. oi M such that supp C 1C for all TV € W and that (x^) is uniformly bounded (in TV). These 
hypotheses will be made throughout this section. 


In this section, we investigate the variational problem which defines and study its minimizer. We 
indeed prove Matytsin’s heuristics iQ outlined in section 6 of |^ . Let us recall the formula obtained in 
(P^ for : 

I^^\nD,d-E) ■=-J 0 {^J■D,^^E) + I^id-E) - inf + 1 - [ x'^dnoix) 


where, for any /r G Mi{IR), 

hid) = ^J x'^dfi{x) - ^S(^). 

J/sidD, ■) is the rate function governing the deviations of the law of the spectral measure of Xm = Dm + Wm 
with a Hermitian (resp. symmetric) Gaussian Wigner matrix Wm and a deterministic diagonal matrix 
Dm = diag((ii,-- - ,dM), (di)i<i<w G , with spectral measure jj,^^ = weakly converging 

towards fin C 'P{1R). It is given (see jl^) by 

Jf}{dD,d) = ^''^x£{Sp^{v,)\v G C([0, 1],V{1R)) : ui = /x}. (2.1) 


if 




- 1 - 00 , 

S°’^{v) := supo< 3 <t<i f ), otherwise. 
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Here, we have set, for any f,g G Cl'^{lR x [0,1]), any s <t G [0,1], and any G C([0, 1],V{1R)), 


= J fix,t)diytix) - J f{x,s)diys{x) 

~ Is I ^ I I 


dxf{x,u) - da;fiy,u) 


x-y 


dvu{x)dvu{y)du, ( 2 . 2 ) 


</,5>L = 


JJ dxf{x, u)dxg{x, u)dvu[x)dv 


(2.3) 


and 




(2.4) 


It can be shown by Riesz’s theorem (see such a derivation in |Q for instance) that any measure-valued path 
V, G C([0,1], Mi(IR)) in {Sf^o < oo} is such that there exists a process k, so that 


1 . 


< / - fc, / - >0,1= 0 


f&cl’\lRxio,i]) 


2. i/Q = go and for any / G Cl'^{Si x [0,1]), any 0 < s < t < 1, 




(2.5) 


Then, it is not hard to show that 


Therefore, Jf 3 {yLz),yt) is given also by 


= ^ >0.1 • 




Jg{yLD,g) = ^inf{< k,k {y,k) satisfies (C)}. 


with (C) the condition 


-L^ {dv±dt) 


( 2 . 6 ) 


(C) : uq = hd, ui = fjL, dxk G C^'^{1R x [0,1])) , and, for any f gC^’ {1R x [0,1]), any s,t G [0,1], 


=< f,k>i. 


The main Theorem of this section states as follows 

Theorem 2.1 Let g,E G ■) < oo} with finite entropy E. Then, the infimum in JpinD, yE) is reaehed 

at a unique probability measures-valued path y* G C([0,1],7^(IR)) such that 

• Mo = Td, Ml = MB- 

• For any t G (0,1), y^ is absolutely continuous with respect to Lebesgue measure ; yl(dx) = pl(x)dx. 
t G [0,1] ^ G V{1R) is continuous and therefore limtj^o Mt = To, limt|i y^ = yE- 
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Let k* be such that the couple {p,*,k*) satisfies (C). Then, if we set 


u; = ^,k: + H^i*iy), 

(p*,u*) satisfies the Euler equation for isentropic flow described by the equations, for t G (0,1), 


dtpU^) = 

9tip*tix)ut{x)) = 


-9x{Pt{x)ul{x)) 


(2.7) 

( 2 . 8 ) 


in the sense of distributions that for all f G x [0,1]), 

J J dtf{t,x)dpt{x)dt + J jdxf{t,x)ut{x)dpt{x)dt = 0 

and, for any e > 0, any f G with fig := {{x,t) G Six [0,1] : pl{x) > e}, 

f (2M*(x)5t/(x,i) + {u^{x)‘^ - n'^plixf) dx:f{x,t)) dxdt = 0. 


(2.9) 


If we assume that {pu,pe) o,re compactly supported probability measures, we additionnally know that [p* ,u*) 


are smooth in the interior of TIq, which guarantees that (2.'i) and (2.8) hold everywhere in the interior of 
LIq. Moreover, TIq is bounded in Six [0,1]. Furthermore, there exists a sequence {(jf)e>o of functions such 
that 




and if we set 


then 


Pf{x):=TT + ^(5a;^‘')^)2 


- d,,^l^^ fdp;{x)dt + Y J {p*{x)-pl{x) f {p*t{x) + pl{x))dxdt < e. 

As a consequence, if we let n((x) = ul{y)dy, which should be thought as the limit in Hi{pf{x)dxdt) of 
the sequence (2“^(/)'^)g>o, we find that it satishes, in the sense of distributions in Oq, 

d,ut = -^id,utr + ^{prr, 

which is Matytsin’s equation [T^ . 

The (non trivial) existence of solutions to the Euler equation for isentropic flow (|2.7|), (|2.8|), is a conse¬ 
quence of our variational study. The uniqueness of the solutions to these equations could be derived, under 
some additional regularity properties, from a convexity property of our rate function 5^^. Even when such 
solutions are not unique, we know that our minimizer is unique due to a convex property of Sfj_^ which is a 
consequence of its representation of Property Ell) below (see Property El)- 
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Property 2.2 Let he G {J^(fj-E, ■) < 00 } having finite entropy S. Then, 

1) For any V & C{[0,l\,V{St)), if{v,k) verifies (C) and ufix) = dxkt{x) + Hvfix), 

S^oir) = \ J(ut{x)fdvt{x)dt + ^J^J{Hvt{x)Ydvt{x)dt - - S(/io)). 

2) Consequently, we can write Jp under the following form 

J 0 {hd, 1 ^e) = J{ut{x)fdfit{x)dt + J^ J{Hfit{x)f dn*{x)dt-{J:{fiE)-^{nD))^ { 2 . 10 ) 


with {fi*,u*) as in Theorem 2.1. Note here that fjf{dx) = pt{x)dx for t S (0,1) and p* G L^{dxdt), so 


that 


J {Hp*t{x))^dptdt= y ^ J {pl{x)fdxdt. 


3) As a consequence. 


Pe) = — 


[ [ {u*t{x))'^dpt{x)dt + [ [ {Hpt{x))'^dpl{x)dt 

Jo J Jo J J 

+ S(^d)) + i J x'^dpoix) + ^ J x'^dpsix)-inf Ip. 


( 2 . 11 ) 


In Q, a similar result was announced (see formulae (1.4) and (2.8) of Q). However, it seems (as far as 
I could understand) that in formulae ( ^.10 2.11 ) of [Q, the first term as the opposite sign. But, in 11, 
formula (2.18), the very same result is stated. 


Let us also notice that the minimizer pt has the following representation in the free probability context. 
Let {A, t) be a non-commutative probability space on which an operator D with distribution pE, 8 x 1 operator 
E with distribution pE and a semi-circular variable S, free with {D,E), live. Then, there exists a joint 
distribution of {D,E) such that {Pt)te[o,i] is the law of a free Brownian bridge 

Xt=tE + {l- t)D + ^t{l - t)S. 

The isentropic Euler equation which governs p* hence partially specify the joint law of {D,E). More 
specifically, for any t G (0,1), our result implies that for any p G IN, 

r Nir{UENU'DM) „ r 

lim / {tUEj^U* + {1- t)DN)P —ET, - dm>l{U) = T{{tE + {1 - t)D)P) = / x^du^ix) 

F^\Dn,Em) j 

if vl is the unique compactly supported probability measure such that, if S' is a semicircular variable free 
with {D, E), for any p G IN, 

j x^dptix) = T [{tE + {1- t)D + ^/t{T^^S)P^ = J xPdEt[^at(i_t){x) 

where 0 denotes the free convolution and as the semicircular variable with covariance 6 . 
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2.1 Study of 

Hereafter and to simplify the notations, = /tq and with some probability measures (/j,o,/j.i) 

on IR. We shall in this section study the rate function and show that it achieves its minimal value on 
{h' G C([0, 1],V{1R)) : vi = ^ 1 } at a unique continuous measure-valued path /r*. 

2.1.1 Sfj,g achieves its minimal value 

Recall that for any probability measure fio G 'P{IR), is a good rate function on C([0, 1],V{1R)) (see Theo¬ 
rem 2.4(1) of 1^). Therefore, the infimum defining J^(/xo, /ri) is, when it is finite, achieved in C([0, 1],V{1R)). 
We shall in the sequel restrict ourselves to (/ro,^i) such that J/3(/ro,/ii) is finite. 


2.1.2 A new formula for S,, 


In this section, we shall give a simple formula of S^g{v) in terms of w. = Hv + dxk, and v when {k,v) 
satisfies (C). We begin with the following preliminary Lemma 

Lemma 2.3 Let (/ro,Mi) ^ {r ^ V{IR) : S(^) > —00} and v. G C([0,1],7^(5?)) such that eq = fj,o, ei = 
and V, G {5^0 < cx)}. Then, for almost all t G (0,1), vt{dx) <C dx and 

3 


I J {Hnt{x)fdEt{x)dt = Yj^ j 


dvtix) 

dx 


dxdt < 00. 


lo J 'J JQ 

The idea of the proof of the lemma is quite simple ; we make, in the definition of , the change of variable 

f{x, t) f{x, t) — f log \x — y\dvt{x). However, because (x, t) —!■ f log \x — yldvtix) is not in C^’^{1R x [0,1]) 

in general, the full proof requires approximations of the path e, and becomes rather technical. This is the 


reason why I defer it to the appendix, section 4.2. We shall now prove the following 


Property 2.4 Let (^O)Mi) & id ^ V[IR) : S(^) > —00} and u, G C([0,1],P(7R)) such that uq = ^oi vx = y,i 
and V, G {S/j-o < o®}- Then, if {v, k) satisfies (C) and if we set Ut := dxkt{x) -I- Hvt{x), we have 

= \ J {ut{x)fdvt{x)dt + \ J {Hvt[x)fdvt{x)dt - ^(S(^i) - E(/ro)). 


Proof. 


Let us recall that (y, k) satisfying condition (C) implies that for any / G C^’ {IR x [0,1]), 


J f{x,t)dnt{x) - J f{x,s)dEs{x) = j J dyf{x,v)diyy{x)ds 

'■ dxf{x,v) - dxf{y,v) 


1 

+ 2 


x-y 


di'v{x)dvx{y)dv 

+ J dxf{x,s)dxk{x,s)diys{x) ( 2 - 12 ) 


with dxk G Lf{dvt{x) x dt). Observe that by |^, p. 170, for any s G [0, 1] such that Vs is absolutely contin¬ 
uous with respect to Lebesgue measure with density ps G Ly,{dx), for any compactly supported measurable 





function d^fi-, s), 


dxf(x,s) - dxf[y,s) 


x-y 


dvs{x)dvs{y) = 2 / dxf{x, s)Hvs{x)dxds. 


Since by Lemma for almost all s G [0,1], Vs{dx) <C dx with a density ps G L^{dx) we conclude that, in 
the sense of distributions on ]R x [0,1], (|2.12|) implies 


^sPs “t” dxiUsPs) — 0: 

i.e for any compactly supported / G C^’°°{]R x [0,1]) vanishing at the boundary of JRx [0,1], 


[ [ {dsf{x,s) +Usdxfix,s))psix)dx = 0. 

Jo Jm 


lo JIR 

Note here that, by dominated convergence theorem, we can equivalently take / G C'^'‘^{]R x [0,1]). 
Moreover, since Hi/ belongs to LJ{dvs x ds) by Lemma 2.3, we can write 


) —< fc,/c >Q i — 


n {us{x))^diys{x)ds + [ [ {Hvs{x)Ydvs{x)d. 
R Jo J M 

nl 


2 / / Hvs{x)us{x)dvs{x)ds 

Jo Jm 


(2.13) 


(2.14) 


We shall now see that the last term in the above right hand side only depends on (/tq, pi). The only difficulty 
in the proof of this point lies in the fact that a;, s G IR x [0,1] ^ / log \x — y\dvs{y) is not in Cl'^{M x [0, Ij). 

However, following Lemma 5.16 in [^, if □ denotes the free convolution (see for a definition), if for 
any i5 > 0, (t^ denotes the semicircular law with covariance (5, and if denotes the field corresponding to 

S(t/i|T]cr^) - S(t/orT]gj) = 2 / f H{iXs\T]js){x)ul{x)di/s\T]js{x)ds. (2-15) 

Jo Jm 

It is well known that 6 S(r|~i~|cr^) is continuous (see ||2^, Theorem 2.1 for the lower semicontinuity and 

use the well known upper semi-continuity). Moreover, if Xg is a random variable with distribution i/g and S 
a semicircular variable, free with Xg, living in a non commutative probability space (M, r), by Theorem 4.2 
in H], the field is given, Vs\T\jS almost surely, by 

uf = T{dxk{Xg, s)|Xs -I- VSS) + Hi/g\T]Ts. 

Consequently, 

{i/g\T\crs){x)Ug{x)di/g\T}Ts{x) = r (dxkjXg, s)H{i/g\T}T 5 ){Xg + +t (H{i/g\T}js){Xg + . 

Moreover, by Voiculescu |^, Proposition 3.5 and Corollary 6.13, if Vg{dx) = pg{x)dx G L^{dx), 

limr (^iHii/g[^s){Xg + VSS) - Hvg{Xg)f'^ = 0. 
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Therefore, for any such s € [0,1], 


lim / H{i'sV+\Js){x)ul{x)di's[T\j 5 {x) = / Hvs{x)us{x)di's{x). (2-16) 

<5^0 JM ^ ^ JR 


Note that by Lemma 2.3, this convergence holds for almost all s G [0,1] since p. G L 3 (dxdt). Finally, by 
Propositions 3.5 and 3.7 of for any s such that Hvs is well defined. 


Hiy,[^siXs + VSS) = T{Hiy,{Xs)\Xs + VSS) 


so that for any (5 > 0 


T ({Hiys\T]js{Xs + VSS))'^'^ < T {{Hi's{Xs) f) . 

Therefore, dominated convergence theorem and ( |2.16| ) imply that 

lim / / Hi'sr+}x 5 (x)u^^(x)dvs[T}Js(x)ds = / Hvs{x)us(x)dvs(x)ds. 

Jm ' Jo Jr 

Thus, (2.15) extends to (5 = 0 which proves, with (|2.14|), Property 2.4. 


2.1.3 Uniqueness of the minimizers of 

We shall use the formula for obtained in the last section to prove that 


Property 2.5 For any {pQ,pi) G V{R) with finite entropy E, there exists a unique measures-valued path 
p* such that 


Jf3{po,h-i) = ^inf{S';,„(y) 


= Ml} = 


In the following, p* shall always denote the minimizer of Property 2.5 and d^k*,u* its associated fields. 
Proof.According to the previous section, the minimizers of 5'^^ also minimize 


S{u,p) = J— [ [ {pt{x))^dxdt + [ [ {ut{x)fpt{x)dxdt 

3 Jo Jr Jo Jr 

under the constraint dtpt + dx{ptUt) = 0 in the sense of distributions, pt > 0 almost surely w.r.t Lebesgue 
measure and f pt(x)dx = 1, and with given initial and terminal data for p given by 


limpt(x)da: = po{dx), 


Irm pt{x)dx = pi(dx) 


where convergence holds in the weak sense (with respect to bounded continuous functions) and is simply 
due to the fact that is finite only on continuous measure-valued paths. 

Let m = up he the corresponding momentum. In the variables {m,p), S{p,m) reads 

S(rn,p) = 'J— [ [ {pt{x))^dxdt + [ [ 

3 Jo Jr Jo Jr pt{x) 

with the convention § = 0, whereas the constraint becomes linear 
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dt{pt{x)) + dx{mt{x)) = 0, pt{x)dx G V{1R) Vt G [0,1], \impt{x)dx 


Poidx), 


lim pt(x)dx 

tn 


pi{dx). 


We now observe that 5 is a strictly convex function. Indeed, if (m^, p^) and {m^,p^) are any two couples 
of measurable functions in {5 < oo}, it is easy to see that for any a G (0,1) 

dl.S{am^ + {1 — a)m'^,ap^ + {1 — a)p^) = 27r^ [ [ {pl{x) — Pf{x))^{apl{x) + {1 — a)pf{x))dxdt 

Jo Jm 


+2 




(pl(x)mUx) - pUx)mUx))^ 

(aplix) + {1 - a)p^{x))^ 


Hence, d^S{am^ + (1 ~ a)m?,ap^ + (1 ~ ^)p^) > 0 for some a G (0,1) unless for almost all t G [0, 1 ] 


i/\ 2 /\ ii/\ (x') 

Pt = Pt (x) = Ptix), and [x) = 1 = 


Kx) 


pi(x) pKx) 


= uf{x) pt{x)dxdt a.s. 


In other words, S is strictly convex. By standard convex analysis, the strict convexity of S results with 
the uniqueness of its minimizers given a linear constraint, and in particular in More precisely, from the 
above, the minimizer p* in Jp is defined uniquely for almost all t G [0,1] (and then everywhere by continuity 
of p*) and its field u*, or equivalently dxk*, is then dehned uniquely p’l[dx)dt almost surely. 

I 


2.2 A priori properties of the minimizer fi* 

In this section, we shall see that the minimizer p* has to be the distribution of a free Brownian bridge when at 
least one of the probability measure po or pi are compactly supported, the other having finite variance (since 
we rely on |^’s results). To simplify the statements, we shall assume throughout this section that both 
probability measures are compactly supported. This property will unable us to obtain a priori properties on 
the laws of the minimizers, such as existence, boundedness, and smoothness of their densities. 

2.2.1 Free Brownian bridge characterization of the minimizer 

Let us state more precisely the theorem obtained in this section. A free Brownian bridge between p^ and pi 
is the law of 

At = (1 - t)Xo + tXi + v^t(l - t)S (2.17) 

with a semicircular variable S', free with Xq and Xi, with law pq and pi respectively. We let FBB(/ro,/ii) C 
C{[0,1],V{1R)) denote the set of such laws (which depend of course not only on po,pi but on the joint 
distribution of (Aq, Ai) too). Then, we shall prove that 

Theorem 2.6 Assume p^, pi compactly supported. Then, 

j0{po,Pi) = = po,i^i = pi} 

= ^inf{S(i^) ; ly e FBB{po,pi)}. 
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Therefore, since FBB{fj,Q, fii) is a closed subset of C{[0,1], V{1R)), the unique minimizer /i* in the above 
infimum belongs to FBB{fio, ^i). 

The proof of Theorem |2.6| is rather technical and goes back through the large random matrices origin of Jjs. 
We therefore defer it to the appendix. 


2.2.2 Properties of the free Brownian motion paths 


As a consequence of Theorem 2.6, we shall prove that 


Corollary 2.7 Assume /rg and compactly supported. Then, 

a) There exists a compact set K C IR so that for all t G [0,1], pf,{K‘^) = 0. For all t G (0,1), the support 
of pf is the closure of its interior. 

b) pf{dx) <C dx for all t G [0,1]. Let pl{x) = 

c) There exists a finite constant C (independent oft) so that, pf almost surely, 

p:{xf + {Hp*{x)f<{t{i-t))-^ 


and 

K(x)|<c(t(i-t))-i 


d) {p*,u*) are analytic in the interior of LI = {x,t G fR x [0,1] : pl{x) > 0}. 

e) At the boundary of Lit = {x G IR : Pti^) > 0}j for x G Lit, 


\pUxrd,,p*t{x)\ < 


1 


47r3t^(l — ty 


P*t (x) < 


4:Tr^t'^{l — ty 


{X - Xo )3 


if xo is the nearest point of x in Lit. 


Consequently, the minimizer p* may only have shocks at the boundary of its support. 

Pro of. This corollary is a direct consequence of Theorem ^ and we shall collect these properties for any 
free brownian bridge law. Indeed, let {A, t) be a non-commutative probability space in which two operators 
Xo,Xi with laws po and pi and a semicircular variable S, free with {Xo,Xi), live. We assume throughout 
that Xo and Xi are bounded by C for the operator norm (i.e po{[—C, Cy) = pi{[—C, Cy) = 0). 

Let pt be the distribution of 


Xt — tXi + (1 — yXo + \/t{l — t^S. 


Clearly, since S is bounded by 2 for the operator norm, Xt is bounded by C + 2 for all t G [0,1]. Thus, 
proposition 4 in Q finishes the proof of a). Following Voiculescu (see Proposition 3.5 and Corollary 3.9 in 
p^), the Hilbert transform of pt is given, /rj-almost surely, by 

Hpt{x) = r((2v^t(l - t)) S\Xt) 


with t( \Xt) the conditionnal expectation with respect to Xt, i.e the orthogonal projection on the sigma 
algebra generated by Xt. We deduce that since S is bounded for the operator norm by 2, ^(-almost surely. 


12 









1 






Further, following Q, the stochastic differential equation satisfied by Xt shows that, for any twice continu¬ 
ously differentiable function / on _/R, 

Mt(/) = Mo(/) + ^ y J J dns{x)dns{y)ds + J j dxf{x)dxks{x)dyLs{x)ds (2.18) 

with k the element of L^{dyLsds) given by 


d^ks{x)=T{ ^^ 

5—1 

Hence, 

Ut := d^kt + H^^t = + T{V^t{l-t)~"s\Xt) = t{X, -Xq+ ^^Z^ S\Xt). (2.19) 

t 1 — t) 

Therefore, //^-almost surely, 

( 2 . 20 ) 

Moreover, by Biane’s results 0, we know that, for t G (0,1), y-t is absolutely continuous with respect to 
Lebesgue measure. We denote by pt its density. Then, we also know that for all t G (0,1), /rt-almost surely, 

pt{xf + {Hptf{x) < (2-21) 

Let us mention the regularity properties that {pt)t^{o,i) will inherite from its free Brownian bridge formula. 
If t't denotes the law of tXi -|- (1 — t)XQ, we have, following Biane |^, corollary 3, that if we set 

v{u,t) = mi{v>0\f -—< (t(i _ 

J [u — xy + 

= inf{u > 0|t (((tXi -I- (1 - t)Xo - uY + v'^y^) < (t(l - t))“^}, 


then 


while 


ip{u, t) = u + t(l ~ J 

Hpt{ip{uy)) = 


{u — x)dvt(x) 

{u — xy + v{u, ty ’ 

(u — x)dvt{x) 


Ptitpiuy)) = 


(u — xy + v{u, ty ’ 

v{u, t) 


Trt{l — t) 

From these formulae, we observe that is analytic in the interior of ft since ip' is bounded below by a 
positive constant there (see Biane ||^, p 713 and the obvious analyticity in the time parameter t G (0,1) ), 
it is clear that p is in Q,. Hence, the weak equation ( 2. IS] ) is verified in the strong sense in H and we find 
that in O, ut{x) = pt(x)~^f^ dtpt{y)dy is C°°. 

At the boundary of fit = {x : {x,t) G fl}, Biane (0, corollary 5) also noticed that 


\ptix)‘^dxdt{x)\ < 


1 


47r3t2(i _ ty 

with Xq the nearest point of the boundary of fit from x. 


Pt{x) < 


- ty 


{x - Xo) 3 
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2.3 The variational problem 

We now turn to the analysis of the variational problem defining Jp ; we shall prove that 

Property 2.8 Assume that /tq and p,i are probability measures on ]R such that S(/ro) and S(/ri) are finite. 
Then, the path p* G C([0, \],V{1K)) minimizing Jp(po,pi) satisfies ; 

1) pI = po and p\= pi. 

2) For any t G (0,1); Tti^x) <C dx. Let {pt)t^(o.i) denote the corresponding density. By continuity of 
p*, pl{dx) = pl(x)dx converges towards po (resp. pi) as t goes to zero (resp. one) in the usual weak sense 
on V{]R). 

3) p* is characterized as the unique continuous measure-valued path such that p^ = po and p* = pi and, 

for any v G {>S'^o < satisfies (C) and vi = pi, we have, with u = d^k + FIv, 

2u’l{utdvt — u)dpl) — j(uf)^(di/t — dp^) J— dpf)]dt > 0 (2.22) 



4)As a consequence, {p*,u*) satisfies the Euler equation for isentropic flow described by the equations, 
forte (0,1), 


dtP*t {x) = 

dt{pt{x)ut{x)) = 


-dx{P*tix)ut{x)) 
-dx{^P*t{x)u*t{xf - yPt(a;)^) 


(2.23) 

(2.24) 


in the sense of distributions that for all f G x [0,1]), 

J J dtfit,x)dpl{x)dt + J J dxf{t,x)ul{x)dpl{x)dt = 0 

and, for any e > 0, any f G with flj := {{x,t) G IR x [0, 1] : pl{x) > e}, 

f {2u^{x)dtf{x,t) + {ul{xY - ■k'^pI{xY) da: fix,t)) dxdt = 0 . 


(2.25) 


Let us now assume that {po,pi) are compactly supported. Then 

5) i 2.2f^ ) is true everywhere in the interior o/flo. Moreover, ( 2.2t ) can be improved by the statement 

that ^ 2 1 

[ [ Utix){dtfit,x)-\-ut{x)da:fit,x))dpt{x)dt=^ [ [ iplix)fda:fit,x)dxdt 

Jo Jm 3 Jo j]R 

for all f G Cl'^ilR x (0,1)). 

6 ) There exists a sequence i(j)'^)e>o of functions such that 


+ >0 


and if we set 


plix):=Tr Adt(j)^-\-4, fida:(j)^) )^ 
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then 


j {ul{x)-d^'^i^fdi4{x)dt + ^ J {pt{x)-pt(x) f{p*t{x)+p'l{x))dxdt<e. 


Discussion 2.9 Matytsin JS] noticed that if we set 

f{x,t) = ul{x) + in p*{x), 


then the Euler equation for isentropic flow implies that f the Burgers equation. Hence, if one assumes that 
f can be smoothly extended to the complex plan, we find by usual characteristic methods that for z € W 


f{f{z,0)t + z,t) = f{z,0) 

and therefore, setting G+{z) = z + f{z, 0) and G-{z) = z — f{z, 1), we see that our problem boils down to 
solve 

G+ o G-{z) = G- o G+{z) = z 

with 3(G'+)(a:) = Trpo(x) and 9(G'_)(a:) = —7rpi(a;) if po and pi are the densities of pQ,p,i respectively. This 
kind of characterization is in fact reminiscent to the description of minimizers provided by P. Zinn Justin 
However, such a result would require more smoothness of (p*,u*) than what we proved here. 


Proof of Property 2.8 : By property 2.4, we want to minimize 

S{p,u) := J J {utix))"^pt{x)dxdt + ^ J J {pt{x))^dxdt 

under the constraint (C’) : 


dtPt + dx{utpt) = 0, \i^pt{x)dx = po, limpt{x)dx = pi 

and when pt{x)dx G V{IR) for all t G [0,1]. To study the variational problem associated with this energy, I 
know essentially three ways. The first is to make a perturbation with respect to the source. This strategy was 
followed by D. Serre in but applies only when we know a priori that {p*,u*p*) are uniformly bounded. 
Since this case corresponds to the case where po,pi are compactly supported, we shall consider it in the 
second part of the proof. One can also use a target type perturbation, which is a standard perturbation on 
the space of probability measure, viewed as a subspace of the vector space of measures. This method gives 
(3) in Property |2.8| as we shall see. The last way is to use convex analysis, following for instance Y. Brenier 
(see Q, section 3.2). We shall also detail these arguments, since it provides some approximation property 
of the field u*, as described in Property |2.S|.6). 

We begin with the target type perturbation. In the following, we denote {p*,u*) the minimizer of S 
under the constraint (C’). Let [p,u) G {S < oo} satisfying the constraint (C’). Then, for any a G [0,1], we 
set, with m = pu and m* = p*u*, 


= (1 — a)p* + ap, m“ := (1 — a){p*u*) + a{pu) := p°‘u°', u°‘ = (m“/p“). 


It is then not hard to check that S{p°‘,u°‘) < oo for all a G [0,1]. Moreover, by the convexity of ^ : a —> 
{pf{x))~^{mf{x) fl+3~^TT^{pf{x) fl for all admissible (p, m), (p', m'), we see that a~^{(j){a) — (j){0)) decreases 
as a ^ 0 showing, by monotone convergence theorem the existence of daS{p°‘, u“)(0+) and 
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9qS'(/o“, m“)(0^) = J [—{u*)^ p* — p + 2mu* +'!t‘^{p*)^{p — p*)]dxdt 

= J [2m* (m — TO*) — {u*Y{p — p*) + 7r^(p*)^(p — p*)]dxdt. 


Hence, for any {p,u) G < oo}, we have 

9Q,S'(/r?“, u“)(0^) = J[2u*{m — m*) — (u*)^{p — p*) + Tr‘^{p*)^{p — p*)]dxdt > 0 (2.26) 


Reciprocally, since S is convex in (p, to) , we know that 

5(p“, m“) > Sip*,u*) + 5a5(p“,M“)(0+)a 

so that ( 2.26| ) implies that S{p°',u°') > S{p*,u*) for all a G [0,1] and (p,u) G {S' < oo}. Hence, (2.26) 
characterizes our unique minimizer, which proves Property |J[3). We can apply this result with 

p = p* + edx4>, TTi = TO* — edt4> 

for some (j) G e > 0, such that dx4>{-T0) = dx4>{-A) = Oj insuring that S{p,u) has finite entropy. 

This yields the second point of Property 2.8.3). Conditions at the boundary of the support can also be 
deduced from ( 2.26 ), but they are hardly understandable, since the conditions over the potentials (j) become 
more stringent. 

To prove the last points of our property which concerns the case where (po, pi) are compactly supported, 
we follow D. Serre |22 and Y. Brenier [^. 

The idea developped in [^ is basically to set at(x) = a{t,x) = (p}(x),p}(a;)uj'(x)) so that div(at(a;)) = 0 
and perturbe a by considering a family 

a® = Jg{a.Vx,th) o g = Jg{p*{dth + u*dxh)) o g 

with a diffeomorphism p of Q = [Oj 1] ^ IR with inverse h = g~^ and Jacobian Jg. Such an approach 
yields the Euler’s equation ( 2.25| ) of Property |2.8| (use the boundedness of {p*,u*) obtained in Corollary ^ 
to apply theorem 2.2 of |Q). Moreover, since we saw in Corollary 2.7.d) that p* and u* are smooth in the 
interior of Hq, ( 2.24 ) results with Property |2.8| .5). 

We now developp convex analysis for our problem following [^. By Corollary 2.7.a), we see that there 
exists a compact K such that /ij (iC^) = 0 for all t G [0, Ij. We set Q = K x [0,1] and E = Cb{Q) x Cb{Q). 
For any continuous functions F,^ G E , we set 


= \F{x,t) + {^^r\idxdt 


if F + (^)^ > 0, on Q, and +oo otherwise. For any (/r, M) G E', let us consider 

a*(/i, M) = sup{ / F{x,t)p{dx,dt) + / ^{x,t)M{dx,dt) — a{F,^)}. 

Jq Jq 


16 

















It is not hard to see that a*{fi, M) < oo iff ^ is non negative, M is absolutely continuous w.r.t /r (to prove 
both cases take = 0 and argue by choosing d) wisely) and /x is absolutely continuous w.r.t Lebesgue 

measure with density in L^{dxdt) (in this case, take d* = 0). Moreover, if we denote ^{dx,dt) = pt{x)dxdt, 
M{dx,dt) = mt{x)dxdt, it is not hard to see that a*{p,M) = S{p,m). Now, let 

P{F,^) = I F{x,t)pf{x)dxdt + I ^{x,t)ul{x)pl{x)dxdt 

JQ JQ 

if there exists (p G such that 

F{x, t) + dtpix, t) = 0, d>(a;, t) + dxPix, t) =0 

for all {x,t) G Q, and is equal to +oo otherwise. We consider 

l3*{p,M) = sup{ / F{x,t)p{dx,dt) + ( ^{x,t)M{dx,dt) — I3{F,^)} 

JQ JQ 

Then, (3* is infinite unless Jq F{x,t){p{x,t) — pl{x))dxdt + fQ^(x,t)(M(x,t} — ml{x))dxdt = 0 for all 
{F, di) G E. Therefore, dtp. + d^M = 0 in the sense of distributions, f p{x,t)dx = 1 for almost all t G [0,1] 
and \imt 10 p{x,t)dx = dpo{x), Yantai p{x,t)dx = dpi(x). As a consequence, 

inf{a*(/x, M) +/3*(/x, M)} = inf{5(p,m) : (p, m) satisfies (C’) and = 0 Vt G [0,1]} 

= 2inf{5'^o(jx) : vi = pi} + {Y,(po) - ^{pi))'■= Z{po, Pi) 


where in the last line we have used Property 2.4 and Corollary 2.7.a). 

Observe that a, (3 are convex functions with values in ] — oo,oo]. Moreover, there is at least one point 
{F, B) G E, namely F = —1, d> = 0 for which a is continuous for the uniform topology on E and (3 finite 
(this is the reason why we need to work on a compact set K instead of IR). Thus, following ||^, by the 
Fenchel-Rockafellar duality theorem (see theoreme 1.11 in [^), we have 


mf{a*{p, M) + P*{p, M), {p,M)GE'} = sup{—a(—F, —d)) —/3(F, d>) : (F, $) G Tl} 

and the infimum is achieved. More precisely, 

Z{po, pi) = snp{J dt4>t{x)plix)dxdt + J dx(t>t{x)rrii{x)dxdt - j (^dtptix) + ^dx4>t{x)'^'^ dxdt} 

where the supremum is taken over p G C^’^(K x [0,1]) such that dt<p+ > 0. As a consequence, there 

exists a sequence of functions in Cl’^{K x [0,1]) such that 


> 0 


and 


Ut{x)'^dp*{x)dt+Y J pl{x)^dxdt< J {dtPHx) + Ut{x)dx(p’') dpt{x)dt-^ j 


>t{x) + -dxptix)^ 


+e^ 
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for all e > 0, which implies, if we set 




J {u^ (x) — dx Ydp^{x)dt < -n^ j pKx)"^p^{x)dxdt — j pl{x)^dxdt — ^ J pl{xYdxdt + 

= j{p*t{x) - pl{x)f {2pl{x) + pl{x))dxdt + (2.27) 


which completes the proof of the Property. 


3 Applications to matrix integrals 

In physics, several matrix integrals have been of interests in the 80’s and 90’s for their applications to 
quantum helds theory as well as string theory. We refer here to the works of M. Mehta, A. Matytsin, A. 
Migdal, V. Kazakov, P. Zinn Justin and B. Eynard for instance. Among these integrals, are often considered 
the following : 

• The random Ising model on random graphs described by the Gibbs measure 

pg,^gidA,dB) = ^_e^tr(AB)-^tr(p,(A))-Artr(p.(B))^^^^ 

^ I sing 

with the partition function 

ZLng = y e'^tr(AB)-Mr(P,(A))-wtr(P.(B))^^^5 

and two polynomial functions Pi,P 2 . The limiting free energy for this model was calculated by M. 
Mehta ^ in the case Pi{x) = P 2 {x) = x'^ + gx'^ and integration holds over Hn- However, the limiting 
spectral measures of A and B under were not considered in that paper. A discussion about this 

problem can be found in P. Zinn Justin ||^ . 

• One can also define the Potts model on random graphs described by the Gibbs measure 

^Potts j—2 

The limiting spectral measures of (Ai, • • • ,Aq) are discussed in ^ when Pi = gx^ — x^ (!). 

• As a straightforward generalization, one can consider matrices coupled by a chain following S. Ghadha, 
G. Mahoux and M. Mehta given by 

^chain i—2 

q can eventually go to inhnity as in |l^. 
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• Finally, we can mention the so-called induced QCD studied in Q. It is described, if A = [—q,q]^ C 


e A) — —/ e 

^QCD igA 


tT{UjAi+eM;Ai) 


2D 


Y\_dm^{Uj)Y\_e dAi 

j—1 iGA 

where iej)i<j< 2 D is a basis of . The description of the limit behaviour of the spectral measures 
of Ai, • • • , Aq is given in [Q in the case q = oo. We impose periodic boundary conditions at the 
boundary of the lattice points A. 


In this section, we shall study the asymptotic behaviour of the free energy of these models as well as describe 
the limit behaviour of the spectral measures of the matrices under the corresponding Gibbs measures. 

The theorem states as follows 


Theorem 3.1 Assume that Pi(x) > CiX^ + di with Ci > 0 and some finite constants di. Hereafter, /3 = 1 
(resp. P = 2) when dA denotes the Lehesgue measure on Sn (resp. Hn)- Then, 




I sing 


Fpotts — 


- + v(Q) - v) - - 2 inf Igiy) 

2 2 ueV{R) 

- d inf 


(3.1) 


(3.2) 


F, 


chain — 


Ar2 ^chain 

N^oo iV^ 


inf 1 / 3 ( 12 ) 

i=i ,;=2 ^ i=i ^eV{lR) 


(3.3) 


Fqcd = 




2D 


-inf{^/i,(P) - --^S(^,)} - 2D|A| inf 1 / 3 ( 12 ) 


ieA 


(3.4) 


Remark 3.2: The above theorem actually extends to polynomial functions going to infinity like x^. However, 
the case of quadratic polynomials is trivial since it boils down to the Gaussian case and therefore the next 


interesting case is quartic polynomial as above. Moreover, Theorem 3.3 fails in the case where P, Q go 
to infinity only like x^. However, all our proofs would extends easily for functions Pfs such that Pfix) > 
a|xP+'^ -I- b with some a > 0 and e > 0. 


Theorem ^ will be proved in the next section, but merely boils down to a Laplace’s (or saddle point) 
method. 
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We shall then study the variational problems for the above energies. We prove the following for the Ising 
model. 


Theorem 3.3 Assume Pi{x) > ax^ + b, P 2 {x) > ax^ + b for some positive constant a. Then 

• The infimum in Fising achieved at a unique couple {p.Ajh-B) of probability measures. 

• {p,A, Pb) are compactly supported measures with finite entropy S. 


Let be the minimizer of 5^^ on {vi = pb} as described in Theorem Then, 

(pA, PB, is the unique minimizer of the strictly convex energy 


C{p,nu,p*,m*) := p{Pi - ^x^) + iy{P 2 - ^x"^) - ji^{p) + ^{v)) 


+ - 


P f f {m*t{x)) ^ 
) J P*t ( 2 :) 


dxdt + — J J pl{x)^dxdt 


Thus, we find that {pA, Pb, p"^^^ are characterized by the property that for any {p, u, p* ,m*) £ 

{£ < 00}, 


^(-Pi - - IJ-a) + j{P 2 - ^x‘^)d{iz - pb) 

'f // ^‘^^\^-v\dLA{v){dp-dpA){,x)-^ J J log\x-y\dpB{y){dp-dpB)ix) 

J - p^^^) + TT‘^{p^^^f{p* - p^^^)]dxdt > 0 


• (p^^^satisfies the Euler equation for isentropic flow with pressure p{p) = —^p^ in the strong 
sense in the interior of LI = {{x,t) G Mx [0,1] : pf^^{x) ^ 0} and satisfy the conclusions of Property 


2 .i. 



• Moreover, 

1 1 j Hi.) - _ 

j {P[{x) -X- {x))h{x)dpA{x) 

and 


j [P^ix) — X — {x))h{x)dpB{x) 

for all h G cl (M). 



Remark 3.4: The last point becomes 

Hpa{x) = P[{x) — X — {x), PA s..s.,F[pb{x) = P^lyX) — X — {x), Pb Si.s. 

as soon as HpA and HpB are in L^[pa) and L^{pb) respectively. 

For the other models, we unfortunately loose obvious convexity, and therefore uniqueness of the mini- 
mizers in general. We can still prove the following 
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Theorem 3.5 • For any given yii, there exists at most one minimizer [^ 2 , ■ ■ • > M 9 ) w Fpotts but unique¬ 

ness of /ii is unclear in general, except in the case q = 3 The critical points in Fpotts are compactly 
supported, with finite entropy E. 

Let {pti, • • • , fiq) be a critical point and for i G {2, ■■■ ,q}, denote ,u^) the unique minimizer described 
in Theorem \2.^ with p^ifix) = pi{dx) and p\{dx) = pfidx). Then 

P'fix) = qx + - ^{q-3)Hpi{x) 

i^2 

in the sense of distributions on supp{fii) and 

P'{x) = X - ^u\{x) - ^Hpfix), 2<i<q 
in the sense of distributions on supp{pi). 


• There exists at most one minimizer in Fchain,- The minimizer {pi,--- , pq) is compactly supported 
with finite entropy E. The critical points {pi, • • • , pq) in Fchain are such that for i G {2, • • • , q\, W It 
is such that if we denote the minimizer described in Theorem ^ with pg^dx) = pi-i(dx) and 

p\(dx) = pi{dx), we have 


P[{x) = X + — ^Hpi{x) and P-{x) =2x — ~ 2 < i < q. 

in the sense of distribution in supp{pi) and supp{pi) respectively. 


Again, uniqueness of the critical points in Fqcd Is unclear in general, except in the case D = 1 where 
uniqueness holds. In this case, the minimizer pi is symmetric, yielding pi = p for all i G A and the 


unique path {p,u) described in Theorem 2.i with boundary data {p,p), satisfies Uq{x) = —u^(x) and 


P'{x) — 2x — fiuQ^x) = 0 pa.s 


3.1 Proof of Theorem 3.1 


The proof of Theorem 3.1 follows a standard Laplace’s method. We shall only detail it in the Ising model 
case, the generalization to the other models being straightforward. 

Let P, Q be two polynomial functions and define, for N G JN, A]\f{P, Q) G IRU {+ 00 } by 

A^(P, Q) = J exp{—A^tr(P(A)) — NtT{Q{B)) + NtT{AB)}dAdB 

where the integration holds over orthogonal (resp. Hermitian) matrices if /3 = 1 (resp. (3 = 2). 

We claim that 

Lemma 3.6 Assume that there exists a,cG , and b,dG IR such that 


P(x) > ax^ + b and Q{x) > cx'^ + d, for all x G M. 
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Then, we have 


sup {-fi{P) - i^iQ) + + -2 inf 

N^oo t,,uGV{m) 2 ^er(lR} 

Remark here that the result could be extend to P{x) > ax^ + b and Q{x) > cx^ + d with ac > 1 but that 
the Gaussian case being uninteresting, we shall use the above and simpler hypothesis. 

Proof. 

Observe that for any e > 0, 

^ |«'~ (r^g) I + 4" ( 1 1 ) 

< Ve ((triv(A‘‘)) " (trjv(S^)) = + (trAr(R'‘)) ^ (trArA^) = ^ 

< ^/e (trAr(yl‘‘) + trAr(B"‘) + trjv(A^) + trAr(B^)) 

Therefore, if we set 

h-^sing{dA,dB) = ^ exp{-A^tr(F(A)) - Ntv{Q{B)) + Ntv{AB)}dAdB 

and 

1 /exp{-iVtr(P(A))-iVtr(Q(i?))+iVtr(^^)}dAdi? 

- ATm) -' 

= log ~ . 

we get 

AAr(e) < ^log/xf^i„g(exp{v^fVtr(A'‘ + A2) + V^fVtr(B^ + R2)}) 

< ^ log/xf,,„g (exp{<jyifVtr(A4 + ^2) + q^,mr{B^ + B^)}) 

where we used Jensen’s inequality with <7 > 1. Now, under our hypothesis, and since 2|yli3| < A"^ + B^, it 
is clear that if q^/e is chosen small enough (e.g smaller than a A c), the above right hand side is bounded 
uniformly. Hence, we take q = 2 a/\c,/e obtain 

limsupAAr(e) < C-s/e (3.5) 
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with a finite constant C. Moreover, for any e > 0, we can use saddle point method (see 0 for a full rigorous 
derivation) and Theorem 1.1 of ||^ to obtain 


lim -^log [ exp{—Ntr(P(A)) — Ntr(Q(B)) — NtY{- 
JV—>oo J ] 


A 


B 


1 49 1 ny)}dAdB 

'l + eA^l + eB^'^ 

with (j)e{x) = (1 + ex^)~^x and /r o (j)~^{f) = fj,(f o cj)^). Thus, (^) results with 


/3 


J™ l^logA 7 v(^’,Q) = linr sup {-/r(P) - i/(Q) +o ,()/, i/o ,()/) +^(E(^) + E(z9))} - 2inf/; 3 . 

Moreover, we can prove as for ( ^ ) that for any ly such that ^i{x'^) < M and I'ix^) < M, 

(/r o ^ o ) - /(«^9 ) I < C(M) v^. 


Using the fact that 


as well as 


E(^) + E(z9) < C{^i{x^) + i>{x^) + 1) 


for some finite constant C, we see that the supremum above is taken at /x, v such that and v{x'^) are 

bounded by some finite constant depending only on P, Q. Hence, we can take the limit e going to zero above 
and conclude. | 


3.2 Proof of Theorem |3.3| and |3.5 

3.2.1 The Ising model 

Let us recall that 

Fising + 2 inf Ii3{v) = - inf{/x(Pi) + v{P2) - v) - ^E(/x) - ^E(zx)} 

v(iV{R) 2 2 


Observe that since v) <2 ^^J^{x'^) + 2 ^v^x"^), the minimizer (/x^, ^b) in the above right hand side is 

such that 

dA{Pi - ^x'^) + /xb(P 2 - - ^E(^ 9 i) - < -Fising - 2 inf Ipi^v) < 00 . 

2 2 2 2 ner{lR) 

Hence, since Pi — and P 2 — 2“^a;^ are bounded below under our hypotheses (for well chosen a), we 

conclude that Y,{^a) and 'S{ij,b) are bounded below and hence finite. Further, if 2ni (resp. 27 x 2 ) is the 
degree of Pi (resp. P 2 ) for ni, 7 x 2 > 2, we also see that 

fiAi.x^"'^) < 00 , ^b{x‘^^^) < 00 . (3.6) 
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Thus, we can use Property to get 




Fising + 2 inf^ Ifiiv) = - inf |/r(Pi - + u{P 2 - - ^(S(fi) + E(i/)) 


f3 


+ T { / 

4 (u*,Ai*)G(C),,,,, Jq 


{J^ j u*t{x)^d^l*t{x)dt + j Hi4{xfdi4{x)dt) 


1 


1 


/3. 


inf <! ^i(Pl -^x) + v{P 2 - 7 :x )- -(S(/r) + E(z/)) 

(u* 


J u*t{xfdfj.*t{x)dt + Y J J p*{xfdxdt 


:= — inf L{fj,,iy, p*,u*) 

(u* ,pi* ) G(C)^,u 

where {u*,fi*) G means that in the sense of distributions 

dtP*t + dxiPtul) = 0, Im p* (dx) = p, lim /i* {dx) = i' 

and we have used in the last line that when the above infimum is finite, is absolutely continuous with 


respect to Lebesgue measure for almost all t G [0,1] and with density p* G L^^dxdt) (see Lemma 2.3). 


Observe that if L{p,v, p*,u*) = C{p,v, p*,m*) with in* = p*u*, £ is a strictly convexe function of 
{p,i/, p*,m*) (recall that —E is convex, see for instance) and that the constraint (C')^_iy is linear in the 
variables {p, p*,m*). Therefore, the above minimum is achieved at a unique point {pA, Pb, P^^^ 

We now perform a measure type perturbation to characterize the infimum. Take (/i, v, p*, m*) G {£ < oo} 
and set, for a G [0,1], 

(/r“, i^°‘, /9“, m“) = a{p, v, p*, m*) + (1 - a){pA, pb,P^^^, 

Then, we find that we must have 


J (Pi - ^x‘^)dip - Pa) + j {P2 - ]^x^)d{v - pb) 

” 2 // ^°s\x - y\dpAiy)idp - dpA){x) - ^ J J log\x-y\dpBiy)idp-dpB){x) 


+ ^ J[2u^^^{m* - - {u^^^f{p* 




Taking p = p^ and b = pB, we see that {p ^^^must satisfy Property Now, if p{dx) = pA{dx) + 

dx(j)o{x)dx, v{dx) = pB{dx) + dx4>i{x)dx and m* = — dttj), Pt = p^^^ + dx4> with (ji G C^’°°{]Rx [0,1]) 

such that , „ 

f + edt(j)f , ^ f {dt(t)f 


. „ „ , -dxdt < oo, / ———dxdt < oo, 

pA-,B^Q P + edx(p JpA^B^Q dx(p 


(3.8) 
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we obtain by (3.7) 


(Pi - ]^x^)dx(j)Q{x)dx + y (-^2 - ^x^)dx(j)i{x)dx 

^og\x - y\d^iA{y)dj;(l)o{x)dx - ^ j J \og\x - y\dfiBiy)dx(l}i{x)da 
^ J[2u^^^dt4> — {u^^^Ydx(j> + Tr‘^{p'^^^)^dx<l)]dxdt > 0 


(3.9) 


which becomes an equality if tp is supported in il. = {{x,t) G JR x [0,1] : ^ 0} by symmetry. If we 

assume that is sufficiently smooth, in particular continuously differentiable with respect to the time 

variable around t = 0 and t = 1, we can use integration by parts to see that 

J[2u^^^dt(l)- {u^^^ fdx4> + n‘^{p^^^ fdx4’]dxdt > 2[Juf^^dx(l)tdx]l 
yielding that there exists two constants li,l 2 such that 

Pi{x) - ^ J log\x - y\dpA{y) -‘2Tlo^^ix) = h pA a.s (3.10) 

P 2 ix)-^x'^J\og\x-y\dpBiy)-‘2Tlf^^{x) = h Pb a.s (3.11) 

Pi{x) - ^x'^ - ^ J log\x - y\dpAiy) - 2 Uq^^{x) > h if a: G supp(/rA)° 

P 2 {x) - ^ J log\x- y\dpB{y) - 2nf^^(x) > h H x e supp(^b)° 


Such a result would generalize the usual equations obtained in the one matrix case. However, since we could 
not prove such a regularity property of {p^^^ we shall now obtain a Schwinger-Dyson type formula 

following iQ, theorem 2.15 and proposition 2.17, to obtain a weak form of ( [l. 10 ),(3.11). Let us briefly recall 
the ideas in the case (3 = 2 (the case /? = 1 being similar), which is based on an infinitesimal change of 
variables. 

If, in we change A ^ A + N~^h{A, B) with some smooth bounded functions h of two non- 

commutative variables (take for instance h belonging to the set CCst{(P) of Stieljes functionals defined in 
[Q, ||](see also its definition in appendice 4.1), it turns out that, due to Kadison-Fuglede determinant formula 
(see |§, the proof of theorem 2.15 and proposition 2.17) 

^ J ^tl-{hiA,B){-P[{A)+B))+N-AmtY{DAh{A,B)}+0{l}-NtT{Pi{A)+P2{B)-AB)j^j^^Q 

with Da the non commutative derivation with respect to A given by 

DA^hg) = DAh xl(g)g + /i(g)lx DaQ, V/i, g G CCst{(P), DaB = 0, DaA = 1 (g) 1. 

Therefore, we can find a finite constant C{Ji) such that for any e > 0 


yLug ® p^^HDaHA, B)) + p^^\{-P[{A) + B)h{A, B))\ < e) > 1 - (3.12) 
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with the empirical distribution of A, B defined by 


= tiN{h{A,B)), \/h e CCst{(P). 


Of course, the same type of formula holds when A is replaced by B. It is not hard to see that is tight 
under for the topology described in p, corresponding to the CC;,t((r)-weak topology (see Q for proof 


Sing 


of similar statements). Let r be a limit point. Taking, for e > 0 and 5 > 0, h{A,B) = (1 + 5A^) Pj{A)(l + 
eB‘^)~^ with j{x) — ni<i< ^{zi — x) ^ for some Zi G([\]R and n G JN, and p large enough {p larger than half 
the degree of P{) so that DAh{A,B) G CCst{^ 0 CCst{(C) and (1 + 5A^)~^{P[{A) — B){\ + eB'^)~^j{A) G 


CCst((I^, we deduce from (3.12) that r must satisfies for any e,S > 0 and p large enough, 


> t(Da(1 + 6A^)-Pj{A) X 1 0 (1 + eB^)-^) = r((P{(A) - B){1 + 6A^)-Pj{A){l + eB^)-!). (3.13) 


Similarly for any e,S > 0, and p large enough, 

T^TiDB{l + 5B^)-Pj{B) X {1 + eB^)-^) = TiiP^iB) - A)il + SA^)-Pj{B)il + eA^)-^). (3.14) 

Now, by (|3.6| ), PHA) — B and P^iB) — A belongs to L^{t) so that we can let 6, e going to zero to conclude 
by dominated convergence theorem that 


T(S)T{DAj{A)) = t{{P[{A) - B)j{A)), T(S)T{DBj{B)) = t{{P^{B) - A)j{B)). (3.15) 

We next show that ( ^.15| ) implies that pA and pB are compactly supported when ni > 2 and n2 > 2, and first 
that all their moments are finite. To this end, take j(a;) = ((1 + ex^)~^x)^ = e“" (l + ie~^{x — ie“^)“^)" {x+ 
^g-i)-n fgj. ^ ^ yielding 

PA (^P[ix) ((1 + ex'^y^x)'^^ = T (jiB\A) ((1 + eA'^y^A)'^^ + t 0 r(i:)^j(^)) (3.16) 

with, since Df can be represented in the tensor product space as Df{x,y) = (x — y)~^(f(x) — f(y)), 

n—1 

T 0 T(DAj(A)) = ^PA (((1 + exy~^x)P) PA (((1 + ex‘^)~^xy~^~P) 

p—0 


n—1 


-e ^ PA (((1 + ex^)-y)P+^) PA (((1 + ex^)-yr-P) . 

p—0 

When n is odd, it is not hard to see that we can find c> 0,dn G ]R such that P'{x)x^ > — dn, so 

that we deduce from (3.16) that 


cpA 


1 2ni — 1+n 


1 + ea;^ 


< dn+ 2n snp PA { {— -y)^ + Pa 

p<n V 1 + ea;^ 


1 + ex^ 


PB{\x\py (3.17) 


where we have used in the last line Holder’s inequality with conjugate exponents p, q. We take q = n~^{2ni — 
1 + n), p = (2ni — l)“^(2ni — 1 + n). Similarly, we obtain for pB, and q = n~^(2n2 — 1 + n), p = 
(2712 — l)“^(2n2 — 1 + n), 


cpB 




Pb 


\nq 


1 + ex"^ 


^^A{\X\PY 


(3.18) 
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Now, we have seen that 


< oo, < oo 


so that ( |3.17 ),( 3.18 ) yields 


= supMA(((l + ea;2)-ia;)2"^-i+") 
£>0 


^ 2 rii-l+n^ _ 


) = sup^B (((1 + ex^) ^x) 

£>0 


2 ^- 1 ^^ 2 ra 2 -l+fe^ 


< oo for 2ni — 1 + n < mf := 2 n 2 { 2 ni — 1) 

< oo for 2 n 2 — 1 + k < mf := 2 n 2 { 2 n 2 — 1) 


and then by induction for 2n\ — 1 + n < rrip := — 1), 2 n 2 — 1 + fc < to® := mp_i{ 2 n 2 — 1) for all 

p >2. Since 2ni — 1 > 1 and 2ni — 1 > 1, rrip and to® go to infinity with p, which proves that pA and ps 
have finite moments of all orders. 


As a consequence, we can extend by dominated convergence theorem ( [1.16 ) to polynomial functions (i.e. 
take e = 0) resulting with 


PA {P[{x)x^) = T {t(B\A)A'-) + Y,pa {x^) pa . (3.19) 

p —0 

and a similar equation for the moments of pB- Let us write P[{x) = + X]p =2 Cipx‘^^^~'P, P^ix) = 

with ai > 0, /3i > 0. Setting a„ = \pa{x'^)\ and 6„ = \pb{x'^)\, we deduce that 


2ni n—1 ^ ^ 

—1+n ^ ^ ^ I^2ni —p+n ^ ^ ^p^n—l—p H” ^qn (3.20) 

p—2 p—0 

2 n 2 n—1 ^ ^ 

/^1^2ni —1+n ^ ^ ^ |/5p|^2ni —p+n ^ ^ ^ ^pbn — l—p H“ bqnClp (3.21) 

p—2 p—0 

with conjuguate exponents (p, q) to be chosen later. 

Now, we make the induction hypothesis that for some R G , for some m G IN, 


Up < RFCp, bp < BFCp, for p < TO 


with Cp the Catalan numbers given by 


p-i 


Cp — CnCp-i-n, Co — 1. 


n—0 

Of course, up to take R big enough, we can always assume that to > 2ni V n 2 - Now, plugging this hypothesis 
into ( 3.20 ),( 3.21 ) with to + 1 = 2ni — 1 + n and q = mn~^, we obtain 

2 ni 

aia2„i-l+n < ^ \ap\R^^^-P+^C2rA-p+n + i?"0„ + i?"+i(C„)^(q^] + i)^ 


p=2 


2 ni 


p =2 
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where we have used that Cm increases with m. Thus, our induction hypothesis is verified as soon as 


2ni 

\ap\R~P + 

p=2 

2712 

\l3p\R-P + + i?2(l-n2) 

p=2 


< ai 


< /3i 


which is clearly the case for R large enough since we asumed ni A 712 > 2. Since m ^ log Cm goes to 4 as m 
goes to infinity, we deduce that 


limsup ^ log + 4, 

m —»^oo ZTTI 


limsup log < i? + 4, 

m —^00 277i 


implying that fiA and are supported into [—7? — 4, i? + 4] for R finite satisfying the above induction 
hypothesis (plus the condition imposed by the first 2ni V 77,2 moments). 


Let us now go back to (^T5) and notice that since the Stieljes functions are dense in Cc(7R) and P{—t{A\B) 
belongs to it can be extended to j G Cl{IR) ; 

j{x)-j{y) 


x-y 


-dy.A{x)dy,A{y) = t{{P'{ x) - T{B\A))j{x)) 


(3.22) 


Since [pA, Ms) are compactly supported, we can use the conclusions of section 
almost surely. 


2 . 2.2 


We see that 


A^B 


,A^B 


= t(B - A\Xt) + (1 - 2t)HyY^{x) 


so that 


,A^B 


= t{B\A) — X + HyA 


at least in the sense of distribution as in (|3.22|) . Thus, by uniqueness of the solutions to the Euler equation 
given the initial and final data {yA, yB) proved in Property 2.5, we conclude that 


in the sense of distribution that 


h{x) - h(y) 
x-y 


HyA{x) = P[{x) - X- Uq^{x) 


dyA{x)dyA{y) = J (T’((x) - x- {x))h{x)dyA{x) 


for all h G Cl{IR). The second equation is derived similarly and one finds that 

2HyB{x) = P 2 {x) - t{A\B){x) = P^ix) - (x - u^^^{x) - HyB{x)) 


resulting with 

HyB{x) = P^ix) — x + Ui^^{x) 

in the sense of distributions. Note also that by Property the fact that {yA, yB) are compactly supported 
implies that {p^^^, satisfies the isentropic Euler equation in the strong sense in Q. 


28 











3.3 q-Potts model 

In this case, we find that 


Fpotts = inf 

i^2 y^V{In) 

= -inf{/ii(Pi - |a;^)+ ^/z,(Pi - y) 

i^2 

4 Uo J Jo J J 

~ 3)S(/^i)} - g inf I^{v) (3.23) 

4 ^ 4 v^v{IR) 


When q > 3, the above functionnal is not anymore clearly convex in fii since S(tti) is concave. Hence, 
the uniqueness of the minimizers is now unclear. Note however that the Euler-Lagrange term may still 
contain sufficient convexity in /ii to insure uniqueness but simply that the above formula does not show it. 
In the case g = 3, the functional is still convex, and strictly convex in the arguments Therefore, 

uniqueness of the minimizers still holds since if (/x, zz, /i*, it*) and ,u*), we would still find that by 

convexity Jl* = /x* and therefore /x = /Xg = )xg = )x, zz = /x* = /x); = zz. The above formula already shows that 
the critical points satisfy iXi{Pi) < 00 and have finite entropy S. We can also obtain the Schwinger-Dyson 
equations in this case and deduce as for the Ising model that the critical points are compactly supported 


and satisfy the equations of Theorem 3.5 


3.4 Chain model 

In this case. 


Fcha 




i=2 

„2 9 


G'P(IR) 


(3.24) 


= -inf{/xi(Pi-—) + ^/x,(Pi-a;2) 


2=2 


P 


z =2 


inf 


1 / 


inf -f/3(j^) 
4 u&v(E.) 


(3.25) 


Here, we still have convexity and strict convexity on the term coming from _ Hence, uniqueness of the 
minimizers hold. Again, we can prove the conclusions of Theorem j3.5| as for the Ising model. 
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3.5 Induced QCD model 


q 2D „ 

Fqcd = - in f{E - EE inf Ip{u) 

i=l ieAj=l ieA v^Viln) 

2 — 1 2 GA 

[ 7 o J Jo J J 

-2£i|A| inf If 3 (iy) 
uev(JR) 

Again, obvious convexity disappears and uniqueness of the minimizers becomes unclear when D > 1. Unique¬ 
ness of the minimizers still holds when D = 1. Then, clearly Hi = n for all i G A and Uq = —u\ at the 
minimizing path with {p*,u*) the solution of the Euler equation with with boundary data {pL^p). p then 
satisfies 


P'{x) — 2x — I3 uq{x) = 0 


in the sense of distributions in supp(^), which corresponds to the result obtained by Matytsin [Q, (4.3)] 
when /? = 2. Actually, since we can prove as for the Ising model that p is compactly supported, it turns out 
that P'{x) — 2x — Puq{x) is in every Lp{dp) and therefore that P'{x) — 2x — (3 uq{x) = 0 almost everywhere 
in the support of p. 


4 Appendice 


4.1 Free Brownian motion description of the minimizers 

Let us return to the probability aspect of the story. In fact, by definition, if 


= XU + R 


N 


jN 


with a Hermitian (if /3 = 2, otherwise symmetric if /3 = 1) matrix X^ with spectral measure p^ and a 
Hermitian (resp. symmetric) Brownian motion , if we denote p^ the spectral measure of X^, then, if 
pU converges towards a compactly supported probability measure pq, for any pi G 'P{1R.), 

limsuplimsuplogP(c?(/t(^,/ri) < (5) = liminf liminflogP(c?(/t(^, ^i) < 6) = —Jp{po, pi). 

S^o N^oo 5^0 N^ao N‘‘ 

Let us now reconsider the above limit and show that the limit must be taken at a free Brownian bridge. 
More precisely, we shall see that, if r denotes the joint law of {Xq^Xi) (the precise sense of which being 


given below) and / j ,’’ the law of the free Brownian bridge (2.17) associated with (Xo, Aii), 
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limsup limsuplogP(fi(/if^,/ii) < i5) < sup limsuplimsuplogP( max , fj-J ) < (5) 

S^O N^oo ^ xoX-i=Mo ^^0 N^oo ^ l<k<n 

v-l 

TO A ^ =M 1 

for any family {ti, ■ ■ ■ ,tn} of times in [0,1]. Therefore, the large deviation estimate obtained in yields 


limsuplimsup^logP(d(/tf,/ri) < S) < o ^ ^ = fJ-i}- 

6 —*-0 N —*-oo ^ 

The lower bound estimate obtained in |l^ therefore guarantees that 

inf{5'(i/),i/o = Mo, vi = /ri} = inf{S'(/r'^), r o = /ro,To Af ^ = ^i}. 


Such kind of result were already obtained in Q and §]. 

Let us now be more precise. We recall that we can define the joint law of the two matrices , X{^ by 
the family 

/l^i(F)=tr^(F(Ao^,Ai^)) 

when F is taken into a natural set T of test functions of two non-commutatives variables and trAr(A) = 
It is common in free probability to consider polynomial test functions. In Q, bounded 
analytic test functions were introduced for self-adjoint non-commutatives variables. T = CCst{(C) is there 
the complex vector space generated by 


f(ai,a 2 )= n 


l< 2 <n 


Zi - a]Xi - a‘fX2 


where (zi)i<i<n belongs to ((f\IR)", 1 < fc < 2 )”^^ to and is the non-commutative product. 

We shall here use the very same set of functions and recall then that the space 


Mo,i = {tGX* : t{I) = 1, t{FF*) > 0, t[FG) = t{GF)) 

is a compact metric space. We denote by D a metric on AJo.i- Let us recall @] that if one considers the 
restriction /i^ = toX^^ of r to functions which only depends on one of the variables A^, k = 1,2, then /i^ is 
a probability measure on IR (in fact the spectral measure of Xk) and that the topology inherited by duality 
from F is the vague topology, i.e. the topology generated by continuous compactly supported functions. 

Since AJo,i is compact, for any e > 0, we can find M G IN, {Tk)i<k<M so that AIo,i C Ui<k<M{T ■ 
D{T,Tk) < e} and therefore 

limsuplimsup^logP((i(/tf,^i) < 5) < max limsup ^ logP(d(/if, ^i) < (5; r^) < e) 


Now, conditionnally to X^, 


or equivalently 


dAf = dH^ 


Af - Af 
1 - t 


dt 
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= tx^ + (1 - t)X^ + {l-t) [\l - sr^dH^. 

Jo 

Let us assume that fiQ i converges towards r G 7Wo,i when N goes to infinity and that Xi,Xq remains 
uniformly bounded for the operator norm. In particular, converges for any t G [0,1] towards 

1^1 = To {tXi + (1 - t)Xo)-^; 

i^;if)=Tif{tX^ + il-t)Xo)) 


for any test function /. Therefore, Voiculescu’s result implies that converges towards the distribution 
fj,l of tXi + (1 — t)XQ + (1 — f) J*(l — s)~^dSs with a free Brownian motion S, free with tXi + (1 — t)Xo. 
We shall now extend this result in our topology and also control the dependence of this convergence with 
respect to the speed of convergence of the distribution of towards r. 

We shall work below with given (X^,X{^) G {d^fig , fj,o) < <5; ^i) < (5; r) < e}. 

Let, for u <t, X^J denote the process 


pu 

= txf + (1 - t)X^ + (1 - s)-^dH^. 

Jo 


Then, one deduces from Ito’s calculus that for any test function / 

/ r\ r-.N 


(/) + -— 


0 


As As X — y 


ds 


(1 - s)- 


+ Mf{u) 


with a martingale M^{u) such that 


E[ sup (Mf (it))2] < 

[0,i] 


WfWl 

iV 2 


Moreover, it is not hard to check that w < t) is tight in C([0, 1],V{]R)) (see the proof of exponential 

tightness of the spectral process of Xg + given in [^). The limit points < t) (when D{pLgi,T) 

goes to zero) satisfy the equation 


Mjv* (/) — if) + 


a-tf 


y-x* < 8 i 


fix)-f'{y). ds 


'12 ■ 


Z Jo " " x-y (l-s )2 

This equation admits a unique solution, as can be proved following the arguments of 
Taking f{x) = and substracting both equations, we find, with 


or |ia, p. 494. 


A^iR)= sup E[|A^.j(e*«")-Mxi(e*«")|], 


\i\<R 


that for u < t 


nu 

A(r(i?) < A^{R) + 4i?2 / Af (i?)ds + 
Jo 


R 

N 

which yields thanks to Gronwall lemma and taken at u = t, since /xj" = 


sup E[|/x^„(e*«"^) - < (^ + sup ~ 


.R 






\i\<R 


\i\<R 


Therefore, if we dehne the distance dp on V{]R) by 

dp{f,,y')= [ |/x(e*«")-/x'(e*«")|e-"«^dC 
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we have proved that there exists a finite constant C such that for all t G [0,1], 


C 

It is not hard to convince ourselves that df is a distance compatible with the weak topology on 'P{1R). 
Observe now that on {d{fif ,fii) < 6,d{fiQ , fio) < <^}j ^ t®’ tight for the usual weak 

topology so that for any e > 0 we can find k > 0 so that for any r and t G [0,1], D{t, /i^i) < e implies 

+ii-t)x^ ^ 

Therefore, for any ti,-- - S [0,1], for any € {d{fi^,fio) < 6;d{fif,fii) < 6\D{fi^-^,T) < e}, 

Chebyshev inequality yields 

P( max dpifixM^d-lJ > < nC{K + 

l<k<n '= iV 

with = ^Xt the distribution of Xt = tXi + (1 — + \/t(l — t)S when the law of (Xq, Xi) is r. Hence 

for any rj, when k (i.e e) is small enough and N large enough, 

E(imax^di.(A^«,/i[J < ?7|Xf) > i. 

Hence 

E(d(Af ,Mi) < < e) < 2F(d(/if,/n) < d; r) < e, max dF(Ax'^ , < d)- 

l<fc<n ife 


We arrive at, for e small enough and any r G AIo.ii 

limsup^logP(d(/lf,/ri) < 6;D{fi^^,T) < e) < lim sup ^ log P ( max < 5). 

AT^oo ’ x^oo l<k<n 

Using the large deviation upper bound for the law of {jd^,t G [0, Ij) from we deduce 


limsup^logP(d(/xf ,/ri) < (5) < 

N^oo 


-min inf 

2 1<P<M fnaxi<fc<„ Frp<i5 


We can now let e going to zero, and then 6 going to zero, and then n going to infinity, to obtain, since S is 
a good rate function, that 


lim sup lim sup log P {d{fi^, fj,i) < S) < 

S —>-0 N —>-oo ^ 


2 


inf 

t-.toXq =^J,Q 
v -1 

TOX^ =M1 


Since it was also proved in |p3 that 


liminfliminf-^logP((i(/t]^,/ri) < d) > — ^ inf Slv) 


we obtain 
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inf SM = inf 

^"0 = ^0 v -1 

ui=fxi t'.toXq =^o 

TOX~^ =fj.l 

Hence, if FBB(/io, /ii) is the set of laws of free Brownian bridges between /tq and /ii, i.e 

FBB(/ro,/ri) = o = fio,ToX^^ = ^i}, 


we have seen that 

inf{5'(z/),r'o = = inf{5'(z/), j/ G FBB(/ro, Mi)}- 


To finish the proof of Theorem 2.6, we need to show that FBB(mo, Mi) is a closed subset of C([0,1], 7^(fR)) 
so that indeed the infimum is reached in FBB(mo)Mi)- 

Observe here that m’" does depend only partially on t since it only depends on {vj G [0,1]}. Noting 
that 

p 

i^lixP) = Y, t^riPr,piX, - Xo, Xo)) 


r—0 


with Pr^p{X, Y) the sum over all the monomial functions with total degree p and degree r in X, we see that 
only depends on the restriction of r to polynomial functions P € S = {Pr,p-, 0 < r < p < cx)}. Of course. 


Mof = {r|5,T G Mo,i,t(X^p + Y^p) < 2C^P,\fp G N} 


is closed for the dual topology generated by the polynomial functions of S. Here C denotes a common 
uniform bound on Xq and Xi, and we have 

FBB(/io,/ii) = G Mo.il = e TWff}. 

We denote, for k G -Mo’f ^^d t G [0,1], i/f G P(JR) the distribution of tXi + (1 — t)Xo when the joint 
distribution of {Xo,Xi) restricted to S is k. Then, m” = We now show that FBB(mo,Mi) is a 

closed set of C([0, 1],P{1R)), which insures, since ^ is a good rate function onC([0, 1\,P{IR)), that the infimum 
is achieved on FBB(mo,Mi)- Indeed, if /x" is a sequence of FBB(mo,Mi) given by + | crt(i-t),f G [0,1]}, 
the weak convergence of pP implies the weak convergence of k". Indeed, for any p G IN, any t G [0,1], 

p^ixP) = ly^-(xP)+Ptipn^'), i<p-^) 


with a polynomial function Pt- Hence, by induction, the convergence of {xP ))(recall that /x" is 
supported by [—C — 2, C + 2] for any n so that weak convergence is equivalent to moments convergence) 
results with the convergence of (a;^)))pg^v> ^'^d again, since (t'r")„g£V i® supported by [—C,C], with 
the weak convergence of towards some probability measure Vf Since this convergence holds for any 
t G [0,1], we can expend the moments in powers of the time variable to conclude that converges towards 
K G . Again by free convolution calculus, this convergence results with the convergence of /x" towards 

/x” G FBB(mo,Mi)- Hence, FBB(mo,Mi) is closed. 
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4.2 Proof of Lemma |2.3| : 


In Hi (see (2.13) and Lemma 2.10) O. Zeitouni and 1 proved that for any path v G C^([0, V\,V{IR)), there 
exists a path such that 

limsupS'°’^(z/''’^) = Sfj,g{iy). 
e.Ai.0 


This path was constructed as follows. Let be the Cauchy law with parameter e and set * /r be 

the convoluted path with the Cauchy law. Moreover, if 0 = < ^2 < ■ ■ ■ < tn = 1 with ti = {i — 1)A, we 


set, for t G [tk,tk+i[, 


= < + 


jt - tk) I 

A ' 




^tk\ 


Let us therefore consider S'°’^(^'^’^). Because we took the convolution with respect to the Cauchy law, the 
Hilbert transform is well defined, and actually a continuously differentiable function with respect to 

the time variable and an analytic function with respect to the space variable. Henceforth, in the supremum 
defining S'°’^(^'^’^), we can actually make the change of function —> f{t,x) — /log|a; — y\di'l’^{y). 

Observing that, with = vt* for f G {0,1}, 


we find that 



y\~ 


^dvt^{y) 


dvl'^{x)dt 


I - s(^^)), 


= I (S(:.D - S(i/^)) + 11' J {Hvl’^fdi^l^^dt 

+ sup { / fidul - [ /odi/o - [ [ dtftdvl'^dt- \< f,f >° 1 ^a} 

/ec^’b[o,i]xlR) J J do J ^ 

> I (S(i.D - E(i.^)) + 11' I 


where in the last line we observed that 


sup { [ fidvl - f fodiy^ - f [ dtftdvl'^dt- \< f,f >°1 ^a} 
/GCj'b[ 0 ,l]xlR) J J do J ^ 

= sup sup {A [ fidvl - X [ fodv^ - X [ [ dtftdvl'^dt - ^ < fj >°’e^A} 

f GCl’^ilo, 1 ] XM) xgM d J Jo d ^ 

^ 1 ( (f fidJ^l - / fpdi^o - fo f dtftdvl'^f \ 

^ /GC^b[ip]xlR) V < /> / ) 

> 0 


Observing that 


[“J 

/ / {Hvl'^fd^il’^dt = A^ / {HvlJ^dvl^ 

d^ d k=od 


35 






converges since t ^ Hv^ and t ^ are continuous for any e C([0,1], 7^(fR)), we arrive at 


limmf ^ (S(^i") - E(z/^)) + ^ J{Hulfdvldt (4.1) 


Notice that for t G {0,1}, 

) = J log \x - y\~^dPe * vt{x)dPe * Vt{y) = \ j log((2; - yf + e^)~^dvt{x)dvt{y). 

Hence, monotone convergence theorem asserts that 

l™S(r't) = S(t't). 

Now, recall that for any p G Tricomi |Q p. 169 asserts that 

YP(a;)^ = ^iHpf{x) - H{p{Hp)){x ), 

so that 

J {Hpf{x)p{x)dx = y y {p{x)fdx. 

Since, for any e> 0, vl is absolutely continuous with respect to Lebesgue measure with density p^ G L^{dx) 
for almost all t G [0,1], we conclude by Fatou’s lemma that that 

+00 > liminf liminf > - (S(i^i) — E(r'o)) + — / I \\Ta\xvi{pl{x))^dxdt. 

efo Aio 2 a Jq J 40 


Finally, it is easy to see that (0) implies that pt{dx) = pt{x)dx for almost all t G (0,1) and then that pl{x) 
converges towards pt almost surely. Hence, we have proved that 


1 2 /* 1 /* 

Sf,o{v.) > - (S(^i) - E(z/o)) + yy J {ptix)fdxdt. 
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